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3 by  the  high  galac 


Two  pulsars  have  been  observed  at  low  radio  frequencies 
using  narrow  bandwidth  to  avoid  pulse  dispersion  problems. 
Pulsar  PSR1133+16  was  observed  at  26.3  MHz  using  the 
University  of  Florida  large  array;  pulsar  PSR2045-16  was 
observed  at  45  MHz  using  the  Maipu  Radiobservatory  large 
array  (Chile) . Long  term  integration  techniques  were  used 

Both  pulsars  were  successfully  detected  and  the  average 
pulse  intensity  and  integrated  pulse  profile  were  obtained. 
Time  resolutions  of  P/38  were  obtained  for  PSR1133+16  and 
P/62  for  PSR2045-16,  where  P is  the  pulse  repetition  period. 
The  two  observed  pulsars  were  selected  because  of  their  low 
dispersion  measure,  long  period  (>1  sec),  high  pulse  inten- 
sity and  their  location  with  respect  to  the  galactic  back- 
ground. 

Several  other  radio  sources  were  also  observed  and  were 
used  to  calibrate  the  effective  area  of  the  26.3  MHz  antenna. 


Chairman 


CHAPTER  1 
INTRODUCTION 


Pulsars  were  discovered  in  1967  by  S.  Jocelyn  Bell  and 
Anthony  Hewish.  The  first  pulsar  paper  was  published  in 
Nature  in  February  1968  by  Hewish  et  al.  (1)  telling  about 
the  discovery  of  pulsar  PSR1919+21  and  giving  some  of  its 
parameters.  Since  then,  about  321  pulsars  have  been  detec- 
ted and  cataloged. 


According  to  currently  accepted  models,  pulsars  are 
rapidly  spinning  neutron  stars  originating  from  supernova 
explosions.  They  are  believed  to  have  strong  magnetic 
fields  on  the  order  of  1012  gauss,  probably  the  strongest 
known  magnetic  fields  in  the  Universe.  The  neutron  star  is 
believed  to  emit  continuously  a narrow  beam  of  radio  waves 


which  oorotates  with  the  star,  producing  the  observed  pulse 
each  time  the  beam  sweeps  past  the  earth.  The  periods  of 
the  pulses  which  have  been  observed  range  from  about  0.03 
to  4.3  seconds.  So  far,  for  only  2 pulsars  has  it  been 
possible  to  find  the  optical  counterpart  of  a pulsar. 

Several  mechanisms  have  been  proposed  to  explain  the 
origin  of  the  radio  emission  and  how  the  rotational  energy 
of  the  neutron  star  is  converted  into  electromagnetic  wave 
pulses;  this  subject  still  remains  one  of  the  least  under- 
stood aspects  of  pulsars.  Some  of  the  important  parameters 


of  a pulsar  that  need  to  be  measured  before  realistic 
models  of  the  radio  emission  mechanism  can  be  developed 
are  the  spectrum  of  the  radiation/  the  pulse  shape,  peak 
intensity,  and  the  time  variations  of  these  parameters.  A 
large  number  of  measurements  have  been  made  in  the  range  of 
frequencies  between  100  and  2,000  MHz,  but  only  a few  below 
100  MHz  (2) . The  low  frequency  part  of  the  spectrum  is 
particularly  interesting  because  most  pulsars  present  a 
turn-over  of  the  spectrum  at  metric  and  decametric  wave- 
lengths, and  also  because  pulse  shape  and  widths  show  some 
peculiarities  in  this  region.  The  spectral  maximum  for 
most  pulsars  is  in  the  range  60  to  250  MHz;  this  peak  is 
nearly  symmetrical  with  respect  to  the  frequency  at  which 
the  maximum  occurs  (2) . Measurements  of  the  low  frequency 
part  of  the  pulsar  spectra  have  been  made  by  Bash  et  al.  at 
38  MHz  (3),  Bruck  and  Ustimenko  at  10,  16.7,  20  and  25 
MHz  (4,  5),  Craft  and  Cornelia  at  40  MHz  (6)  and  more 
recently  by  Izvekova  et  al.  at  61  and  102.5  MHz  {7,  8). 

Most  pulsars  measured  at  low  frequencies  are  in  that 
part  of  the  sky  which  is  observable  from  northern  hemis- 
phere observatories.  Observations  at  the  lowest  frequen- 
cies have  been  carried  out  in  the  southern  hemisphere  only 
by  Australian  observatories,  at  80  and  150  MHz  (9,  10) . 
Detection  of  pulsars  at  low  radio  frequencies  is  difficult 
because  of  the  deterioration  of  the  signal-to-noise  ratio 
caused  by  the  high  galactic  background  temperatures,  and 
the  low  frequency  spectral  turn-over  exhibited  by  most 


pulsars.  Also  pulse  dispersion  in  the  interstellar  medium 
becomes  more  severe  the  lower  the  frequency.  At  higher 
frequencies,  compensation  can  be  made  for  dispersion,  per- 
mitting the  use  of  relatively  high  bandwidths.  However, 
at  the  lower  frequencies,  it  is  generally  necessary  to 
restrict  the  bandwidth  because  of  dispersion,  further 
reducing  the  sensitivity  of  the  detection  system.  Due  to 
the  intrinsic  short  pulse  durations  of  most  pulsars  (S  150 

another  limitation  on  the  sensitivity  of  the  system. 

Nevertheless,  the  sensitivity  of  the  system  can  be 
improved  by  averaging  many  independent  records  of  the  same 
source.  If  N records  are  averaged,  the  signal-to-noise 
ratio  (S/N)  is  improved  by  a factor  of  /n  . For  pulsars, 
this  can  be  accomplished  by  averaging  groups  of  consecutive 
periods.  Digital  sampling  is  generally  used  in  such  aver- 
aging. Due  to  the  relatively  short  pulsar  periods,  the 
signal  needs  to  be  sampled  at  a fast  rate  in  order  to  have 
enough  time  resolution  to  meet  the  requirements  of  the 
sampling  theorem.  Since  it  is  not  feasible  to  manipulate 
the  large  amounts  of  data  generated  in  the  digitizing  proc- 
ess by  hand,  digital  computing  techniques  must  be  used, 
providing  an  efficient  way  to  average,  process,  store  and 
plot  that  data. 


In  the  past  few  years,  the  development  of  small  and 
inexpensive  microprocessors  has  made  possible  improved 
methods  of  pulsar  data  handling  and  processing.  If  the 


period  of  a pulsar  is  known  to  be  a good  accuracy,  it  is 
possible  to  average  many  periods  (N  typically  in  the  range 
of  a thousand  periods),  greatly  improving  the  S/N  ratio. 

The  University  of  Florida  Radio  Observatory,  located 
in  Dixie  County  about  50  miles  west  of  Gainesville,  and  the 
Maipd  Radioastronomical  Observatory  of  the  University  of 
Chile,  located  20  miles  southwest  of  Santiago,  have  long 
been  engaged  in  a cooperative  program  of  low  frequency  radio 
astronomy.  Although  most  of  the  effort  over  the  years  has 
gone  into  the  study  of  Jupiter's  radio  emission,  other 
decameter-wavelength  sources  are  also  investigated.  Each 
of  the  two  observatories  has  a very  large  antenna  array,  in 
addition  to  many  smaller  radio  telescope  antennas.  The 
Florida  array  consists  of  640  dipoles  operating  at  a center 
frequency  of  26.3  MHz,  while  the  Chile  array  has  528  dipoles 
and  operates  at  a center  frequency  of  45  MHz.  Both  arrays 
were  believed  to  be  suitable  for  observing  pulsars,  although 
previous  attempts  to  use  them  for  this  purpose  had  not  been 
successful. 

The  goal  of  the  research  on  which  this  thesis  is  based 
was  to  develop  a low  frequency  pulsar  data  acquisition  and 


reduction  system,  and  tc 
important  parameters  of 
two  large  arrays  (i.e.  t 


one  pulsar  with  e 


CHAPTER  2 

DATA  ACQUISITION  AND  REDUCTION 

Two  pulsars,  PSR1133+16  and  PSR2045-16,  were  selected 
(from  the  lists  in  Appendix  A)  for  observation  because  of 
their  low  dispersion  measure,  relative  long  period 
(>1  second),  high  pulse  intensity  and  favorable  location 
with  respect  to  the  galactic  background.  Pulsar  PSR1133+16 
was  observed  at  26.3  MHz  using  a bandwidth  of  B kHz,  and 
pulsar  PSR2045-16  was  observed  at  45  MHz  with  a bandwidth 

The  detected  IF  signals  of  the  pulsar  were  recorded  on 
one  channel  of  a magnetic  tape,  and  a time  code  generator 
signal  on  the  other  channel.  The  tapes  were  played  back 
(making  a 4:1  slow  down),  the  pulsar  signal  sampled,  and 
the  resulting  values  stored  in  the  NERDC  Amdahl  470  computer. 
The  time  code  generator  signal  provides  the  timing  pulses 
for  the  sampling  and  also  indicates  Universal  Time  (U.T.). 

The  sampling  of  the  signal  was  accomplished  by  a 12  bit 
analog/digital  converter  (A/D)  and  a KIM  I microprocessor. 
Pulsar  pulse  averaging  is  carried  out  by  a computer  program 
which  averages  consecutive  periods  of  the  sampled  pulsar 
signal,  previously  stored  in  the  Amdahl  computer. 


2.1  Data  acquisition  of  Pulsar  PSR1133+16 

The  pulsar  PSR1133+16  was  observed  for  20  nights 
between  April  25,  1979,  and  June  8,  1979,  using  the  Uni- 
versity of  Florida  26.3  MHz  large  array,  located  at  the 
University  of  Florida  Radio  Observatory  in  Dixie  County. 

The  26.3  MHz  array  is  a nearly  filled  rectangular  array  of 
640  half-wave  dipoles  that  covers  about  30,000  m2.  The 
half  power  beamwidth  (HPBW)  is  2.5°  north-to-south  by  6.0° 
east-to-west.  It  is  a phase  steered  array  with  full  N-S 
directional  control  of  the  beam  in  quarter-beamwidth  incre- 
ments, and  a selection  of  8 E-W  beam  directions  spaced  one 
beamwidth  apart.  Phasing  is  accomplished  by  means  of  a 
network  of  Butler  matrices,  hybrid  rings  and  plug-in  phasing 
cables . 

North-south  phasing  of  the  array  for  a source  having  a 
given  declination  is  accomplished  in  two  ways,  a coarse 
adjustment  in  large  increments  of  declination  and  a fine 
adjustment  providing  small  steps.  A change  in  the  coarse 
phasing  requires  the  replacement  of  many  plug-in  cables 
located  throughout  the  array,  and  requires  about  2 or  3 man- 
hours of  labor.  Fine  phasing  for  the  desired  declination  is 
quickly  accomplished  by  means  of  switches  at  a single  loca- 
tion. The  beam  pointing  can  be  altered  in  small  increments 
to  any  declination  within  a strip  roughly  10“  wide  by  means 
of  the  fine  phasing  controls.  Positioning  the  beam  to  a new 
declination  outside  this  strip  requires  a coarse  phasing 
readjustment.  The  right  ascension  of  the  beam  is  of  course 


controlled  by  the  rotation  of  the  earth.  Scans  by  any  or 
all  of  the  8 selectable  east-west  beams  (designated  4E,  3E, 

2E,  IE,  1W,  2W , 3W,  4W)  are  made  by  starting  each  scan  at 
the  appropriate  sideral  time.  A computer  program  called 
ARRAY  developed  by  Dr.  M.  D.  Desch,  is  used  to  obtain  infor- 
mation on  cable  lengths  required  for  coarse  phasing,  switch 
positions  for  fine  phasing,  and  beam  transit  time  across  the 
desired  source  (11) . A more  detailed  description  of  the 
antenna  and  its  adjustment  has  been  given  by  Desch, (12)  and 

In  order  to  use  the  maximum  gain  of  the  array,  only  the 
central  E-W  beam  positions  2E,  IE,  1W  and  2W  were  used,  and 
within  each  beam  only  the  central  portion  of  it  was  recorded 

corded  each  night.  Pulsar  PSRU33+16  was  selected  because 
of  its  relatively  high  pulse  energy  (**-0.75x10”^*’  Wm'^Hz-*  at 
61  MHz),  long  period  (1.188  s) , and  low  dispersion  (DM  - 5 
parsec  cm  ) . It  was  also  easy  to  reach  this  pulsar  using 
the  phasing  the  antenna  had  at  that  time  (phased  to  make 
Jupiter  observations) , without  having  to  completely  rephase 
the  array.  Only  the  fine  phasing  needed  readjustment.  During 
the  period  in  which  the  pulsar  was  observed,  it  was  transiting 
at  about  midnight. 

Only  8 nights  (out  of  the  20  observed)  gave  reliable 
data,  free  from  static  and  stations.  The  end  of  the  observa- 
tion period  coincided  with  the  beginning  of  the  thunderstorm 


season,  at  which  time  it  was  necessary  to  shut  the  array 
down  for  the  summer. 

A block  diagram  for  the  interconnections  of  the  equip- 
ment used  to  observe  the  pulsar  is  shown  in  Figure  2-1.  The 
equipment  basically  included  a Collins  R 390  receiver,  a 

quency  standard,  and  a Systron  Donner  8154  time  code  gener- 
ator/reader. 

The  time  code  generator  (TCG)  was  synchronized  against 
WWV;  this  was  accomplished  by  triggering  an  oscilloscope  with 
the  1 pulse  per  second  (1  pps)  output  of  the  TCG  and  measur- 
ing the  delay  of  WWV  1 PPS.  The  1 PPS  of  the  TCG  was  brought 
as  close  as  possible  to  the  1 PPS  of  WWV  (usually  between 
5 to  20  ms)  by  making  short  interruptions  of  the  1 MHz  signal 
from  the  Tracor  frequency  standard.  The  delay  between  the 

i0.5  ms,  and  from  that  information  the  U.T.  was  obtained, 
except  for  the  propagation  delay  between  Boulder,  Colorado, 
and  Dixie  County  Radio  Observatory.  No  corrections  have 
been  made  for  this  delay  because  no  attempt  was  made  to 
relate  the  epoch  of  this  observation  with  any  other  observa- 

The  measured  drift  of  the  Tracor  time  standard  is  about 
80  us/h;  this  causes  a drift  in  the  epoch  of  about  100  us 
between  pulsar  transits  by  beams  2E  and  2W  (•vih16m) this 
drift  is  small  enough  to  be  noticed,  considering  the  time 


r PSRU33+16. 


resolution  of  '>>31.5 


period)  obtained  for 


The  detected  IF  output  of  the  receiver  was  recorded 
on  channel  1 of  the  Magnecord  at  15  IPS  and  the  recording 
level  for  the  galactic  background  was  set  to  about  7 VU  on 
the  recording  level  monitor,  to  prevent  any  saturation  during 
signal  peaks. 

The  time  code  generator  was  recorded  on  channel  2 of  the 
Magnecord.  This  signal  serves  a double  purpose:  it  provides 

a time  reference  for  the  pulsar  signal,  and  later  on,  in  the 
data  reduction  process,  provides  the  timing  signal  for  the 
KIM  I microprocessor. 

As  a back-up  to  obtain  the  correct  time,  at  the  begin- 
ning of  each  observing  session  the  WWV  signals  were  recorded 
on  channel  1 (the  TCG  also  on  channel  2)  for  about  2.5 
minutes.  The  correct  time  can  be  obtained  by  measuring  the 
delay  between  the  1 PPS  of  WWV  and  the  1 PPS  of  the  TCG. 

The  receiver  was  tuned  to  a center  frequency  of  26.3  MHz 
and  the  bandwidth  set  to  8 kHz.  Using  data  for  higher  fre- 
quencies an  extrapolated  pulse  width  (at  the  base)  of  about 
80  ms  was  obtained.  Pulse  broadening  due  to  dispersion 
(DM  = 5 parsec  cm-3)  is  about  18  ms  (for  a bandwidth 
Av  = 8 kHz) , which  was  considered  acceptable. 

hi Data  Acquisition  of  Pulsar  PSR2045-16 

Pulsar  PSR2045-16  was  observed  for  11  nights  between 
May  15  and  May  30,  1979,  using 


the  Maipfi  Radio  Astronomical 
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Observatory  45  MHz  large  array,  located  about  20  miles 
south-west  of  Santiago,  Chile.  The  45  MHz  array  is  a filled 
rectangular  array  of  528  full-wave  dipoles  with  an  effective 

pendently  phased  subsections  in  the  N-S  direction  and  can 
be  oriented  to  a zenith  angle  of  about  i45°  in  the  meridian 
plane.  The  HPBW  is  2.1*  north-to-south  by  4.6°  east-to-west. 

A multiple  beam  system  is  generated  by  a Butler  matrix.  A 
more  detailed  description  of  the  array  has  been  given  by 
Reyes  (14)  and  May  et  al.  (15). 

Pulsar  PSR2045-16  was  recorded  for  20  minutes  each  night, 
which  is  the  time  it  takes  to  drift  between  the  E-w  3 dB 
points  of  the  beam.  This  pulsar  was  selected  because  of  its 
relatively  high  pulse  intensity  (*vl. 0xl0~2®Wm~^Hz~*  at  34  MHz), 
long  period  (1.96  ms)  and  its  low  dispersion  measure  (DM  = 12 
parsec  cm  ) . A zenith  angle  of  only  18°  at  transit  was  also 
another  favorable  factor  considered. 

The  equipment  used  for  the  data  acquisition  of  this 
pulsar  consisted  of  a modified  General  Dynamics  telemetry 
receiver,  a Magnecord  1022  magnetic  tape  recorder,  a Systron 
Donner  8220  TCG  and  a Rhode  and  Schwartz  quartz  time  standard. 
The  interconnections  are  similar  to  those  shown  in  Figure  2-1. 
During  the  observing  sessions  the  pulsar  was  transiting  at 
dawn,  and  no  interference  was  noticed.  Nine  of  the  eleven 
nights  gave  reliable  data;  the  other  two  nights  were  discarded 
because  of  equipment  malfunctions. 
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A bandwidth  of  10  kHz  was  used;  this  is  the  minimum 
bandwidth  of  the  receiver.  Using  data  obtained  at  higher 
frequencies/  an  extrapolated  pulsar  pulse  width  at  the 
base  of  about  130  ms  was  obtained  at  45  MHz.  The  pulse 
broadening  due  to  the  10  kHz  bandwidth  is  11  ms.  This  was 
considered  acceptable. 

The  quartz  time  standard  is  checked  daily  against  the 
Chilean  National  Observatory  cesium  time  standard;  the  epoch 
is  believed  to  be  known  with  an  accuracy  of  about  ±10  ps. 

The  drift  in  time  is  about  3.6  us/h  (or  about  0.00190  ps/P, 
where  P = 1.96  s,  the  pulsar  period).  The  drift  is  1.2  ps 
in  the  20  minutes  between  the  east-to-west  HPBW.  For  all 
practical  purposes  considered  in  this  research,  this  is 
negligible. 


2.3  Data  Reduction 

Data  reduction  for  both  pulsars  was  performed  using  the 
facilities  of  the  Department  of  Astronomy  of  the  University  of 
Florida  in  GaineBville.  The  basic  configuration  and  inter- 
connections of  the  equipment  for  data  reduction  are  shown  in 
Figure  2-2.  The  schematics  for  the  most  relevant  blocks  are 
presented  in  Appendix  C. 

As  stated  earlier,  the  IF  signal  in  the  receiver  is 
band  limited  to  8 or  10  kHz  and  appropriately  smoothed.  The 
resulting  audio  frequency  signal  is  recorded  in  analog  form 
on  channel  1 of  the  Magnecord  tape  reaorder  at  a tape  speed 
of  15  IPS.  The  time  code  generator  is  recorded  on  channel  2. 
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In  order  to  achieve  the  desired  time  resolution  and  to  be 
able  to  sample  the  data  at  an  appropriate  rate,  the  magnetic 
tapes  were  played  back  at  3 3/4  IPS,  providing  a slow  down 
of  4:1. 

Pulsar  signals  from  channel  1 of  the  Magnecord  are  sent 
through  a diode  detector  and  smoothing  filter  to  the  12  bits 
A/D  converter.  The  filter  time  constant  was  adjusted  to  45 
ms  for  both  the  rising  and  falling  edges  of  the  noise  pulses. 

The  TCG  signal  from  channel  2 of  the  Magnecord  is  con- 
nected to  the  input  of  the  TCG  (working  now  as  a reader)  to 
provide  both  a visual  display  of  the  time  and  a 1 PPS  pulse 
(taken  from  a modified  output).  The  latter  gives  the  syn- 
chronizing signal  to  start  the  sampling  of  the  data  at  a 
precisely  known  time.  It  is  particularly  important  to  know 
the  time  at  which  the  sampling  is  started  in  order  to  be  able 
to  combine  in  phase  2 or  more  blocks  of  data. 

TCG  signal  from  channel  2 of  the  Magnecord  is  connected  also 
to  a phase  locked  loop  (PLL)  circuit.  The  PLL  circuit  locks 
at  a particular  frequency  of  the  signal,  filtering  out  the 
rest  of  the  frequencies  present  in  the  complex  wave  form  of 
the  TCG  signal.  The  1000  Hz  signal  of  the  TCG  is  transformed 
to  250  Hz  due  to  the  4:1  slow  down.  The  PLL  was  carefully 
tuned  to  250  Hz,  and  the  output  then  divided  by  2.  The  resul- 
tant 125  Hz  constitutes  the  timing  signal  for  the  KIM  I. 

A gated  DC  output  from  the  PLL  circuit  was  monitored  on 


recorder.  This  provides  a good  means 
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to  detect  any  loss  of  synchronism  or  dropout  of  the  TCG 
playback  signal.  If  there  is  any  loss  of  synchronism 
between  the  input  and  PLL  signals,  then  the  gated  DC  output 
falls  to  zero,  giving  a visual  indication  to  stop  the 
sampling  process.  If  synchronization  were  lost,  the 
pulsar  signals  would  be  averaged  out  of  phase,  wiping  out 
any  pulse  at  the  end  of  the  integration  process.  A means 
was  provided  for  starting  the  sampling  at  the  next  1 PPS 
pulse  from  the  TCG  occurring  after  a switch  was  thrown. 

This  made  it  possible  to  know  the  exact  time  at  which  sampling 
started.  This  was  accomplished  by  connecting  the  start  switch 
to  the  D input  and  the  1-PPS  signal  to  the  clock  (c)  input  of 
a D-type  flip-flop.  The  output  of  the  flip  flop  gates  the 
125  Hz  signal  in  a NAND  gate  (see  schematic  in  Appendix  C) . 
When  the  start  switch  is  turned  on,  the  next  1 PPS  of  the 
TCG  enables  the  125  Hz  signal  to  go  into  the  KIM  I and  the 
sampling  is  started. 

The  sampled  values  (coded  in  hexadecimal)  are  sent 
through  the  Hazeltine  terminal  to  the  Amdahl  470  computer 
and  stored  there  for  further  processing  and  analysis. 

Since  the  maximum  number  of  cards  that  can  be  stored 
under  a new  data  file  name  is  493,  it  was  necessary  to 
break  up  the  data  into  several  "blocks"  to  meet  this  require- 
ment. The  8 minutes  of  data  obtained  in  one  beam  for  pulsar 
PSR1133+16  were  broken  up  in  2 blocks  of  4 minutes  each, 
generating  about  312  cards  for  each  block  (for  a TSAMP=  02 
and  timing  signal  frequency  = 125  Hz;  see  Appendix  B) . The 


20  minutes  of  data  obtained  each  night  for  pulsar  PSR2045-16 
was  broken  up  into  4 blocks  of  5 minutes  each,  generating 
about  390  cards  for  each  block. 

The  method  of  breaking  up  the  data  into  blocks  proved  to 
be  useful  also  because  one  can  inspect  the  results  of  the 
average  of  each  block  separately,  making  it  possible  to  search 
by  eye  for  a possible  pulse  that  might  appear  at  the  same  loca- 
tion in  the  different  blocks.  It  is  highly  unlikely  that 
interference  transients  would  appear  in  the  same  location  in 
each  of  several  blocks. 

Monitoring  of  the  DC  gated  output  of  the  PLL  circuit  is 
very  useful  since  synchronism  is  sometimes  lost  when  a tape 
is  played  back.  If  that  happens,  then  that  portion  of  the 
tape  is  skipped  and  a new  block  started. 

The  KIM  I program  which  is  necessary  to  digitize  the 
data  is  loaded  into  the  microprocessor  from  a cassette  tape 


2.4  Block  Designation  for  Each  Pulsar 
Many  blocks  of  data  were  generated  in  the  data  reduction 
process.  In  order  to  avoid  any  confusion  with  the  data  and 
to  be  able  to  distinguish  quickly  between  blocks,  a designa- 
tion scheme  was  developed  for  referring  to  a specific  data 
block.  The  designation  format  contains  in  abbreviated  form 
the  pertinent  information  for  that  block.  For  pulsar 
PSR1133+16,  the  following  block  designation  was  adopted: 


— Number  that  designates  the  night 
when  the  pulsar  was  recorded 

— Letter  that  designates  the  block 

— Beam  in  which  the  pulsar  was  recorded 


'Designates  pulsar  PSR1133+16 

m designations  are  as  follows: 

st  beam  east  of  transit 


In  order  to  distinguish  between  the  different  blocks 
within  a given  beam,  a capital  letter  was  used,  usually  A a 
B.  In  case  it  became  necessary  to  break  up  the  data  into 
more  than  two  blocks  (i.e.  due  to  loss  of  synchronism  of  th 
TCG  signal)  C was  also  used. 

The  indicator  for  the  date  of  the  night  on  which  the 
pulsar  was  observed  was  an  integer,  as  follows: 

1 April  25,  1979 

2 April  27,  1979 

2 April  28,  1979 

4 April  30,  1979 


1979 


The  designation  scheme  for  the  blocks  of  data  fr 
pulsar  PSR2045-16  is  slightly  different,  as  follows: 


L„ 


t designates  the  night  wh 
the  pulsar  was  recorded 
— Letter  that  designates  the  block 
within  that  night 
— Designates  pulsar  PSR2045-16 
The  block  designation  is  usually  one  of  the  letters  A 


The  reduced  dat 
called  UF.B0030601.S 

tape  called  PULSAR. 


May  17 
May  18 
May  19 
May  22 
May  23 


May  30,  197 
was  stored  c 


XX,  where  XX 
n recorded  a 


a digital  magnetic 


2.5  Computer  Programs  for  Data  Reduction 
In  order  to  analyze  the  data  and  to  plot  the  values,  a 
set  of  two  computer  programs  called  PULSAR  and  MCAPGM  are 
used.  The  basic  parameters  necessary  to  analyze  the  data  must 
be  given  to  the  PULSAR  program;  they  are: 

N - Number  of  bits  digitized  (=12) 

TPP  = Timing  pulse  period  (ms) 

NPA  = Number  of  pulse  intervals  averaged 
P = Pulsar  period  (s) 

MAG  = Plot  magnification  factor 

The  MCAPGM  program  (using  the  PULSAR  program  parameters) 
takes  the  first  group  of  pulse  intervals  to  be  averaged  (NPA) 
from  the  data  of  a particular  block  and  calculates  the  over- 
all average  and  the  standard  deviation.  These  two  quantities, 
and  the  average  for  each  bin,  are  plotted.  The  program  next 
takes  the  second  set  of  intervals  from  the  data  and  does  the 
same  as  before.  Also  the  first  set  of  intervals  is  averaged 
with  the  second  set,  computing  a new  average  (cumulative)  for 
each  bin,  an  overall  average,  and  a new  standard  deviation 
and  plots  the  values  for  the  bins.  This  process  continues 
including  each  time  the  next  group  of  intervals  (NPA)  . 

The  first  set  of  programs  is  run  once  only  to  check  the 
quality  of  the  data, to  detect  any  problem  and  to  obtain  the  total 
overall  average  for  a block  which  will  be  used  in  a second  run. 

The  second  set  of  programs  called  PULS1T0  and  MCAPGM1 
is  similar  to  the  already  described  PULSAR  and  MCAPGM,  but 


this  time,  the  data  is  taken  from  disk  (where  it  has  been 
previously  stored) , produces  a magnified  plot  and  punches 
cards  with  the  final  values  for  each  bin,  the  total  block 
average  and  the  final  standard  deviation.  This  set  of  cards 
is  used  later  on,  to  stack  2 or  more  blocks. 

Usually  after  the  process  already  described,  the  digi- 
tized data  is  recorded  on  the  PULSAR  magnetic  tape  and 
released  from  the  data  files. 


In  order  to  stack  2 or  more  blocks,  the 
of  each  block  (except  the  first,  which  is  use 
is  shifted  the  proper  number  of  bins  so  that 
averaged  in  phase.  The  number  of  bins  to  be 
block  is  given  by  the  following  relationship: 


id  as  reference) 


respect  to  block  1 
= Time  start  sampling  b 
= Time  start  sampling  b 
= Pulsar  period 
= Time  resolution  (=  jjjj 
= Fractional  part 


a block  i needs  t 


be  shifted  with 


The  fractional  par 
the  period  between  t 
time  of  beginning  o 


of  - 


2 — corresponds  to  the  fraction 
of  beginning  of  block  i and 
irst  block.  The  relationship 


s evaluated  using  an  HP-29C  calculator. 


A last  program  called  P2045  (or  P1133I  reads  the  cards 


CHAPTER  3 
DATA  ANALYSIS 

3.1  PSR1133+16  Data  Analysis 
From  the  8 nights  of  data  reduced  for  pulsar 
PSR1133+16  (26.3  MHz,  Florida),  pulses  were  found  with 
no  ambiguity  on  2 nights.  Pulses  were  found  on  April  28 
and  May  25,  1979,  over  averages  of  340  and  390  periods 
respectively,-  the  results  are  plotted  in  figures  3-1  through 

On  April  28,  after  averaging  over  340  periods  (i-6ln44s) 
a pulse  appears  in  bin  33  with  an  amplitude  of  4,3  o above 
the  mean  value  (Figure  3-1) . The  2 data  blocks  P112EA3  and 
P112EB3  used  to  obtain  this  final  block  have  been  plotted 
separately  in  figures  3-2  and  3-3.  It  can  be  seen  that  a 
pulse  appears  in  both  blocks  at  the  same  location  (bin  33) . 
Even  though  the  pulse  appearing  in  each  block  would  not 
alone  be  considered  to  be  statistically  significant,  the 
fact  that  they  appear  in  the  same  location  for  two  sepa- 
rated sets  of  data  provides  an  excellent  test  of  the  validity 

On  May  25,  after  averaging  over  390  periods,  a promi- 
nent pulse  was  obtained  in  bin  25  with  an  amplitude  of 
3.4  o;  this  pulse  is  shown  in  figure  3-4.  The  same  test  for 
validity  as  before  was  made.  Blocks  P112WA6  and  P112WB6  are 
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plotted  separately  in  figure  3-5  and  3-6.  A pulse  appears 
in  bin  25  in  block  P112HA6,  with  an  intensity  2.6  a above 
the  mean  value;  also  in  block  P112HB6,  the  pulse  is  weaker, 
being  only  2.3  a,  but  is  at  the  same  location  (bin  25). 

Once  again  we  have  the  fact  that  2 pulses  appear  in  the 
same  location  in  two  independent  sets  of  data. 

Prom  a precise  knowledge  of  the  pulsar  period  and  the 
elapsed  time  between  the  April  28  and  May  25  appearances  of 
the  pulse,  it  should  be  possible  to  preserve  its  phase  so 
that  it  would  reappear  in  the  same  bin  on  the  second  date 
as  it  occupied  on  the  first.  However,  this  could  not  be 
done  because  of  the  period  uncertainty.  This  problem  will 
be  addressed  later  when  analyzing  PSR2045-16  data. 

The  pulsar  was  detected  in  only  one  beam  on  each  of  the 
two  nights,  disappearing  when  the  beam  was  switched  to  the 
next  position.  This  phenomenon  will  be  discussed  later  in 
the  conclusions,  when  polarization  effects  are  considered. 

Before  other  tests  made  on  this  pulse  are  discussed,  it 
is  interesting  to  note  that  in  bins  9 and  21  of  block 
P112WA6  (figure  3-5)  two  more  small  pulses  can  be  distin- 
guished. Although  neither  of  them  alone  would  be  considered 
statistically  significant,  they  both  seem  to  appear  again  in 
block  P112WB6 (figure  3-6)  in  the  same  position.  The  pulse 
in  bin  21  leads  the  main  pulse  in  bin  25  by  125  msec 
(resolution  -v3i  msec/bin)  or  approximately  38°  of  phase. 

The  pulse  in  bin  9 is  separated  by  496  msec  from  the 
main  pulse,  which  is  about  150°  (or  0.42  P,  P = 1.88023  sec). 
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The  pulse  in  bin  21  suggests  the  existence  of  a 
subpulse  leading  the  main  pulse.  Several  authors  have 
indeed  reported  for  this  pulsar  a complex  main  pulse  com- 
posed of  2 subpulses , the  leading  subpulse  decreasing  in 
intensity  with  respect  to  the  trailing  subpulse  as  the 
frequency  decreases  (6,  7,  8).  The  separation  between  the 
two  subpulses  has  been  measured  at  higher  frequencies  and 
extrapolation  from  that  data  suggests  a much  closer  separa- 
tion for  low  frequencies  than  the  separation  obtained  at 

and  Ustimenko  (4)  reported  some  observations  made  at  25  MHz 
with  a linearly  polarized  antenna;  they  found  that  this 
pulsar  shows  a variety  of  pulse  shapes  with  rapid  transi- 
tions from  one  to  another.  In  one  plot,  they  show  a pulse 
composed  of  2 subpulses  and  separated  by  '■148  msec.  Despite 
the  poor  resolution  used  (P/16  or  t74  msec) , the  2 subpulses 
are  clearly  defined. 

The  separation  of  the  subpulses  as  a function  of  fre- 
quency has  been  plotted  in  figure  3-8,  including  the  result 
by  Bruck  and  Ustimenko  and  our  result  (i.e.  for  the  bin  21 
and  25  subpulses) . The  plot  suggests  that  either  the  pulse 
separation  increases  very  rapidly  toward  lower  frequencies, 
or  the  subpulse  detected  at  26.3  and  25  MHz  is  a new  one 
which  appears  only  at  low  frequencies.  In  figure  3-9  the 
average  pulse  profile  obtained  in  this  thesis  is  compared 
with  that  obtained  by  Bruck  and  Ustimenko.  The  phase  of  the 


pulse  profile  at  26.3  MHz  has  been-  shifted,  to  align  the 
pulses  in  bin  21  and  25  with  the  pulses  at  25  MHz. 

Although  the  subpulse  in  bin  21  is  much  weaker  than 
that  in  bin  25  for  the  390-period  average,  the  two  are 
almost  the  same  strength  in  the  140-period  average  shown  In 
figure  3-10  (block  P112EA6) . In  the  latter  case,  the 
heights  are  2.68  o and  2.62  o for  the  bin  25  and  21  pulses, 
respectively.  A possible  implication  is  that  the  bin  21 
pulse  varies  in  intensity  much  more  than  does  that  of  bin  25. 

A small  pulse  appears  in  bin  9,  when  blocks  P112WA6  and 
P112WB6  are  averaged  together.  This  pulse  is  also  present 
in  each  block  average  separately  (figure  3-5  and  3-6)  , sug- 
gesting the  possible  existence  of  an  interpulse  leading  the 
main  pulse  by  0.42  P (or  150").  In  Bruck  and  Ustimenko  (4) 
a broad  and  poorly  defined  pulse  seems  to  appear  at  0.5  P 
from  the  main  pulse,  in  the  same  plot  where  the  double  main 
pulse  appears (figure  3-9).  So  far,  no  interpulse  has  been 
reported  for  this  pulsar. 


3.1. a Peak  Flux  Density. 

Density  and  Energy  for  Pulsa 

Several  continuous  radio  source 
to  calibrate  the  array,  to  obtain  th 
flux  density,  and  then  to  obtain  an 
flux  density  for  pulsar  PSR1133+16. 


: PSR1133+16 


es  were  observed  in  ord 
he  minimum  detectable 
estimate  of  the  peak 


n parameters  are 


presented  in  Table  3-] 


Table  3-1.  Radio 


used  for  calibration  purposes. 


Radio  Source 


Declination 


Flux  Antenna 

Density  Beam 

(Jy)  Used 


381  IE 

267  1W 


Positions  and  flux  densities  were  taken  from  the  Clark 
Lake  Survey  at  26.3  MHz  (16).  Of  these  four  sources,  only 
3C409  and  3C433  were  used  in  the  calibrations;  the  3C310  and 
3C315  measurements  were  discarded  because  their  difference  in 
R.A.  is  only  8.8m  (corresponding  to  2.2“).  Since  the  HPBW  of 
the  radiotelescope  is  6.0“  east-to-west , these  two  sources 
cannot  be  resolved. 

The  parameters  of  the  radiotelescope  for  observing  the 
calibration  sources  were: 

Bandwidth  (Av)  =250  kHz 

Frequency  (v)  = 26.3  MHz 

The  minimum  detectable  flux  density  Smin,  deduced  from 
calibrations  using  these  sources,  is  presented  in  Table  3-2. 


Table  3-2. 


Radio  Source 


Galactic  background 


3C409 

3C433 


The  galactic  background  temperature  for  pulsar 
PSR1133+16  is: 

Tbackgr.  = 18<514»K 

All  the  galactic  background  temperatures  given  are  actual 
brightness  temperatures,  corrected  for  losses  in  the  antenna. 

A Hewlett  Packard  noise  generator  was  used  for  the  cali- 
brations, and  a pair  of  type  5722  noise  diodes  served  as  the 
standard  of  comparison  for  the  noise  generator.  Details  of 
the  calibrations  are  given  in  Appendix  E. 

The  minimum  detectable  flux  density  has  the  following 
expression  (17): 


where  n = number  of  averaged  records 

Av  - pre-detection  bandwidth  of  the  rece; 
i = post-detection  time  constant 
Tsyst  = system  temperature 
In  this  case,  the  temperature  of  the  system  i 
dominated  by  the  galactic  background  temperature. 


(3-1),  it 
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A.  In  order  to  be  consistent  it  has  been  considered 
that  the  value  Sm^n  for  both  the  continuous  source  and  the 
pulsar  is  that  which  would  cause  a deflection  of  3 a (i.e. , 
3 standard  deviations)  from  the  mean  value  of  the  galactic 
noise  background,  'l'he  3 o criterion  was  used  to  determine 
whether  or  not  a pulsar  pulse  was  present  at  the  end  of  an 
integration  run. 

B.  It  is  assumed  that  the  sample  rate  is  fast  enough 
compared  with  the  post-detection  time  constant  that  very 
little  information  is  lost  in  this  process;  then  K can  be 
considered  tl.  (Details  about  sampling  rates  can  be  found 
in  Appendix  B . ) 


C.  It  is  assumed  that  the  pulsar  signal  is  unpolarized 
as  in  the  calibration  source.  This  is  probably  not  correct; 
a strong  linearly  polarized  component  has  been  reported  (18) 
The  effect  of  variable  Faraday  rotation  in  the  terrestrial 
ionosphere  would  then  be  to  cause  the  plane  of  the  linearly 
polarized  component  to  be  oriented  (at  different  times)  at 
any  angles  from  parallel  to  perpendicular  with  respect  to 
the  antenna  dipoles.  As  a result,  if  the  polarization  is 
actually  100*  linear,  the  value  of  Smin  as  given  by  (3-2) 
may  range  anywhere  from  too  low  by  a factor  of  two  to  too 
high  by  a factor  of  infinity.  Thus,  no  signal  at  all  would 
be  detected  if  the  linear  polarization  is  perpendicular  to 
the  dipoles,  while  according  to  (3-2)  S . would  not  be 


(3-2) 
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By  making  these  assumptions,  once  the  pulsar  reaches 
a value  slightly  higher  than  3 o at  the  end  of  an  averaging 
process,  the  peak  flux  density  can  then  be  considered  equal 
to  Smin  •oalculated  for  the  corresponding  number  of  periods 
averaged).  Therefore, 

speak  p “ Smin  p = 150  J* 

Pulse  energy  is  usually  defined  as: 

Ep=J  S(t)  dt  (3-3) 

where  P = pulsar  period 

S(t)  = flux  density  as  a function  of  time 
For  pulsars  with  no  interpulse,  the  limits  of  the  integral 
are  restricted  to  the  duration  of  the  pulse  (base  time  dura- 


in  order  to  make  the  computations  simpler,  it  is  usually 
assumed  that  the  pulse  has  a relatively  simple  geometrical 
shape.  Izvekova  et  al.  (7)  assumed  a triangular  shape  for 
most  pulsars  (including  PSR1133+16).  By  looking  at  the  pulse 
shape  obtained  on  May  25,  1979,  one  can  see  that  this  is  a 
reasonably  good  approximation,  considering  only  the  pulse 
centered  at  bin  25. 

Using  this  approximation,  the  pulse  energy  in  equation 
(3-3)  can  be  written  as: 


where  Wfaase  = pulse  width  a 


, corrected  for  dispersion 


receiver 


constant. 


Substituting 


i following  v 


defined  as: 

(3-51 

he  value  for  the 


The  mean  flux  density  of  the  pulse  i 

Using  the  value  obtained  for  E and 
period  P = 1.188  sec,  the  following  value  for  the  mean  flux 
density  is  obtained: 

Smean  = 4 4 45  Jy 

It  is  interesting, now,  to  compare  these  values  with 
values  obtained  by  other  authors. 

In  1973,  Bruck  and  Ustimenko  ( 4 ) observed  the  same 
pulsar  at  frequencies  of  16.7,  20  and  25  MHz  using  a linearly 
polarized  broadband  array  and  reported  peak  flux  densities  in 
the  range  of  10  to  40  Jy. 

In  1968,  Drake  and  Craft  (19)  observed  the  same  pulsar, 
using  the  Arecibo  1000-foot  radio  telescope  at  frequencies 
of  111.5  and  195  Mhz  and  reported  peak  flux  densities  of 
150  Jy  at  111.5  Mllz.  Since  these  observations  were  made  at 
a much  higher  frequency,  they  cannot  be  directly  compared 
with  the  one  at  26.3  MHz,  nevertheless,  the  value  obtained 
at  the  two  frequencies  gives  an  idea  of  the  order  of  magnitude 
for  the  peak  flux  densities.  A value  of  38  Jy  can  be  deduced 
from  the  data  at  61  MHz  given  by  Izvekova  et  al.  (7). 


Regarding  the  pulse  energy  and  mean  flux  density,  more 
recently,  in  1979,  Izvekova  et  al.  (7) , have  given  some 
results  at  61  MHz.  They  obtained  a value  for  the  pulse 
energy  of  0.89x10  Jm  ^Hz  ^ and  a mean  flux  density  of 
0.75x10”^®  Jy,  using  a linearly  polarized  antenna. 

In  1971,  Slee  and  Hill  (9),  using  a circularly  polarized 
antenna,  reported  at  80  MHz  a value  for  the  pulse  energy  of 
2.2x10  ^ Jm  ~Hz  * for  one  night  (with  a median  over  6 nights 


It  is  not  easy  to  compare  the  values  obtained  in  this 
thesis  with  values  obtained  by  other  authors,  because  of  the 
polarization  used  and  the  different  averaging  times  (not 
always  given) . Nevertheless,  the  value  of  pulse  energy 
Ep  = 5.3x10  JnT^Hz  ^ obtained  at  26.3  MHz,  looks  at  least 
comparable  with  the  value  of  E = 2.2x10”^®  Jm-^Hz  ^ obtained 
at  80  MHz  (which  also  corresponds  to  a peak  flux  density  of 


■vllO  Jy) . 

According  to  Izvekova  et  al.  (7),  PSR1133+16  has  a 
maximum  in  its  average  spectrum  at  vmax  = 100±40  MHz. 
Xuz'minetal.  (20),  making  simultaneous  observations  at 
several  frequencies,  found  that  the  instantaneous  spectrum 
changes  from  day  to  day.  For  several  days  it  doesn't  show  a 
maximum  in  its  spectrum  but  rather  the  intensity  seems  to 

This  can  be  considered  as  part  of  the  explanation  of 
why  the  peak  flux  density,  the  mean  flux  density  and  the 
pulse  energy  obtained  at  26.3  MHz  are  higher  than  the  values 
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obtained  at  higher  frequencies:  the  pulsar  seems  to  show 

abnormally  strong  pulses  at  low  frequencies  during  some 
particular  nights,  as  it  may  have  done  in  our  case. 

Another  interesting  aspect  of  this  has  been  pointed 
out  by  Ekers  and  Moffet  (21).  They  observed  PSR1133+16  at 
a wave  length  of  23  cm  (''-1300  Mhz) . From  a histogram  of 
energy  and  number  of  pulses  with  a given  energy,  they  noticed 
a long  tail  in  the  distribution,  towards  the  high  energy  end. 
The  ratio  of  the  mean  value  (0.020xl<T26  jm_2Hz_1)  to  the 
highest  energy  pulse  (*v0.090xl0~2^  Jm~^Hz-2)  is  about  4.5. 

If  one  takes  the  value  for  the  mean  energy  of  5.3*10-2®Hz~* 
at  26.3  MHz  as  an  average  over  a relative  short  period 
(•'■390  PI,  and  compares  it  with  the  calculated  mean  energy  of 
1.6x10  2®  Jm  ^Hz  ^ , obtained  from  the  observations  of  Bruck 
and  Ustimenko  at  2S  MHz  (4)  a value  of  3.3  is  obtained  for 
the  ratio.  Such  a ratio  is  not  improbable  according  to  the 
distribution  function  of  Ekers  and  Moffet.  Pulse  energies 
obtained  reported  in  this  thesis  and  by  several  other  authors 
are  plotted  as  a function  of  frequency  in  Figure  3-11. 

3. l.b  PSR1133+16  Pulse  Width 

Several  factors  can  distort  the  actual  pulse  shape  as 
observed:  the  factors  that  need  to  be  considered  are: 

A.  Effect  of  the  dispersion  over  the  finite  bandwidth 

B.  Effect  of  the  post  detection  time  constant  of  the 
receiver. 


C.  Polarization  of  the  pulses  and  antenna. 

The  first  two  factors  will  make  the  pulse  wider  than  the 
actual  pulse,  but  their  effects  can  be  corrected  for  if  the 
pre-detection  bandwidth,  post-detection  time  constant  of  the 
telescope  and  the  dispersion  measure  for  the  pulsar  are 
known;  a simple  form  needs  also  to  be  assumed  for  the  pulse. 

sider  the  basic  formula  for  the  time  of  propagation  through 
a plasma  with  a dispersion  measure  DM  at  a frequency  v. 

The  time  of  propagation  is  given  by: 


where  k is  a constant  depending  on  the  plasma. 

Differentiating  this  equation  with  respect  to  v,  taking 
finite  At  and  Av  (instead  of  dt  and  dv) , and  expressing  the 
parameters  in  consistent  units,  the  following  expression  can 
be  obtained: 

At  = 8.3*10-3DMAv(v„)-3  ms  (3-7) 

where  DM  = dispersion  measure  (parsec  cm-3) 

Av  = bandwidth  (kHz) 

vn  = center  frequency  (hundreds  of  MHz) 

This  relationship  gives  the  time  it  takes  for  a 
dispersed  pulse  to  sweep  through  a bandwidth  Av  centered  at 
frequency  vB . 

If  a triangular  shape  for  the  pulse  is  assumed,  (3-7) 
is  transformed  into  (3-8)  for  the  half  intensity  level: 


The  correction 


constant  of 


(for  the  half  intensity  level  and  assuming  a triangular  pulse 
shape)  is  given  by  (8): 

To  ■=  0.693t  (3-91 

where  t is  the  post  detection  time  constant  of  the  receiver. 

No  correction  can  be  made  due  to  the  effect  of  polariza- 
tion but  some  discussion  about  its  effect  will  be  given  in 
the  conclusions. 

Consequently,  the  corrected  half  intensity  width  of 
the  pulse  is: 

W0.5  = W0.5  ob  ' 4t0.5  " Tc  (3-10) 

where  WQ  is  the  observed  half  intensity  width. 

Substituting  in  equations  (3-8)  and  (3-9)  the  values 


V.  = 26.3  MHz 
and  T = 11.25  ms, 

the  following  values  are  obtained  for  the  corrections: 


Table  3-3  gives  the  values  for  the  observed  half 
intensity  width  and  the  half  intensity  corrected  width  for 
the  pulses  obtained  on  April  28  and  May  25,  1979.  The  pulse 
width  at  the  base  is  the  width  assuming  a triangular  shape, 
and  is  twice  the  corrected  half  intensity  width. 


47 


3-3. 


pulses  are  indeed  d 


it  is  almost  a certainty  t 

May  25,  1979,  the  date  on  which  the 
obtained,  was  chosen  to  make  the  test.  [ 
using  periods  longer  and  shorter  than  the 
by  ±40,  ±80,  ±160, ±320,  and  ±640  us. 

Periods  shorter  or  longer  than  the  assumed  period  P 
by  ±80  us  should  cause  a 


i ta  were  analyzed 
assumed  period 


s shown  in  Figure  : 


plotted 


390  period  average. 

The  result  of  this  test 
the  relative  averaged  pulse 
as  a function  of  the  assumed  period  length,  taking  as  the 
base  the  period P = 1.188023  s given  by  the  DOPVEL  program. 
It  can  be  seen  that  the  pulse  nicely  peaks  with  the  assumed 
pulsar  period,  and  decreases  sharply  for  shorter  and  longer 
periods.  The  lack  of  symmetry  of  the  plot  might  be  due  to 
the  unsymetrical  pulse  shape. 

The  computer  plots  of  the  results  for  the  different 
periods  used  in  this  test  are  presented  in  Figures  3-14  to 


The  pulse  intensity  as  a function  of  time  has  been 
plotted  in  figure  3-25  for  the  pulse  obtained  on  May  25, 
1979.  This  pulse  was  chosen  because  it  was  the  better  of 
the  two  and  is  well  defined.  This  plot  is  useful  not  only 
to  visualize  the  temporal  behavior  of  the  intensity  of  the 
pulse  but  also  to  check  that  no  interference  or  isolated 
transient  pulse  is  mistakenly  taken  as  a pulsar  pulse. 
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regarded 


double 


previously 


described. 

Each  point  on  the  graph  represents  the  intensity  with 
respect  to  its  corresponding  average  value  of  10  consecutive 
periods,  for  bin  25  (peak  of  the  pulse  at  the  end  of  the 
averaging  process) . The  intensity  scale  is  in  arbitrary 
units  and  the  time  scale  is  in  periods  P given  by  the  DOPVEL 
program  (P  = 1.188023  s).  The  plot  shows  that  the  intensity 
in  bin  25  is  almost  always  positive  with  respect  to  its 
corresponding  10-periods  average  value , and  no  clear  periodic 
variations  appear  over  the  390  periods  plotted. 

3,2  PSR2045-16  Data  analysis 

Of  the  10  nights  of  data  reduced  for  pulsar  PSR2045-16, 
pulses  were  found  on  4 nights:  May  15,  16,  18  and  19,  1979. 

A summary  of  the  dates,  averaged  pulse  deflection  in 
terms  of  standard  deviations  e,  and  number  of  periods 
averaged  is  presented  in  Table  3-4. 

Table  3-4.  summary  of  the  date,  pulse  deflection  and 
number  of  periods  averaged  each  night  for 
pulsar  PSR2045-16 


Pulse  deflection 


Periods  averaged 


May  15,  1979 


540 


415 
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Plots  of  the  data  for  these  dates  are  presented  in 
Figures  3-26  to  3-30.  Figure  3-31  corresponds  to  a plot 
of  a typical  day  when  only  noise  was  present  and  no 
pulse  appeared. 

In  general  the  pulses  were  of  different  shapes, 
except  for  the  last  two  days,  for  which  the  shapes  were 
similar.  Some  have  steep  leading  and  trailing  edges  (as 
on  Hay  16,  Figure  3-28);  others  have  only  a steep  trail- 
ing edge  and  a not  so  well  defined  leading  edge  (as  on 
May  15,  Figure  3-26).  On  May  18  and  19  (Figures  3-29 
and  3-30)  both  pulses  show  a steep  leading  edge  and  much 
less  steep  trailing  edge. 

On  May  16  the  pulse  was  clearly  present  with  an 
intensity  of  3.48  a above  the  mean  value  over  an  average 
of  only  180  periods  (Figure  3-27).  On  the  other  3 days, 
pulses  were  not  seen  clearly  until  the  end  of  the  averag- 
ing process  for  that  night.  For  partial  averages  in 
these  cases,  some  indication  of  a pulse  can  be  seen  at 
the  location  where  the  pulse  appears  after  completing  the 
average,  but  its  intensity  never  rises  significantly 
above  the  noise  level.  No  indication  of  the  existence  of 
a complex  pulse  or  interpulse  can  be  deduced  from  the 
plotted  data. 

Pulsar  PSR2045-16  presents  at  higher  frequencies  a 
complex  pulse,  composed  of  3 subpulses,  with  the  trailing 
subpulse  decreasing  dramatically  in  intensity  with 
decreasing  frequency  (221;  at  610  MHz  the  trailing  pulse 
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rapidly,  and  at  151  MHz  is  about  2/3  of  the  central 
subpulse.  The  intensity  of  the  central  and  leading 
subpulses  seems  to  remain  almost  constant  with  fre- 
quency. No  indication  of  resolved  subpulses  can  be 
obtained  from  the  data  at  45  MHz.  Nevertheless,  the 
difference  in  steepness  of  the  leading  and  trailing 
edges  might  suggest  the  existence  of  an  unresolved 
much  less  intense  subpulse  accompanying  the  main  sub- 
pulse. This  aspect  will  be  discussed  later  when  dis- 
cussing the  pulse  width  and  the  effect  of  polarization. 

3.2. a Peak  Flux  Density,  Mean  Flux 


In  order  to  estimate  the  peak  flux  density  of  this 
pulsar,  several  radio  sources  were  observed  for  calibration 
purposes.  Radio  source  PKS2034-17,  one  of  the  closest 
to  the  pulsar,  was  clearly  recorded  despite  the  steepness 
of  the  galactic  background  (23) . This  radio  source  has 
apparently  not  been  observed  before  at  any  frequency  lower 
than  80  MHz.  However,  an  extrapolated  flux  density  at 
45  MHz  was  obtained  from  the  spectral  index  derived  from 
flux  density  data  at  higher  frequencies,  as  reported  in  the 
literature. 

Spectral  index  can  be  defined  (17)  by  the  relation: 


Density 


3-16 


(3-11) 


where  S = flux  density 


v = frequency 
n - spectral  index 

For  any  two  frequencies  on  the  linear  part  of  the 
relationship  between  S and  i , 


(3-12) 


The  best  available  values  of  the  flux  density  for  this 
source  (24,  25)  are  at  frequencies  of  160  and  BO  MHz: 
si6o  =  S * *  8 * *-8  Jy 
S80  “ 15 

Substitutinq  these  values  into  equation  (3-12) , a value 


of  n = 0.76  is  obtained. 

Assuming  a linear  relationship 
down  to  115  MHz,  the  flux  density 


for  S as  a function  of 
at  this  frequency  is 


S45  = 23  J* 

Applying  the  previously  used  criteria,  the  following 

value  for  the  minimum  detectable  flux  density  at  45  MHz 

was  obtained  for  the  region  near  source  PKS2045-17: 

Smin  = 6 Jy 

The  radiotelescope  parameters  for  observing  the  radio 


source  PKS2045-) 


79 

where  Av  = bandwidth,  t = post  detection  time  constant 
and  n - number  of  records  averaged. 

It  can  be  shown  that  the  brightness  temperature  T 
and  the  frequency  v are  related  in  a similar  way  as  S and 
V are  related  in  equation  (3-12) : 


Using  the  galactic  background  temperature  at  38  MHz 
given  by  Blythe  (26) : 


where  Tp  = galactic  background  temperature  near  the  pulsar 
7*r  = galactic  background  temperature  near  the  radio 
source  PKS20 34-17 

and  the  value  for  the  galactic  background  spectral 
index  given  by  Findlay  (27) , valid  for  the  range  20  to  400 
MHz: 


a galactic  background  temperatures  T 


Using  equation  (3-2)  and  the  appropriate  parameters 
of  the  radio  telescope,  the  minimum  detectable  flux 
density  can  be  obtained  for  the  pulsar  observations. 


telescope  parameters  when  observing 


The  parameters  when  observing  the  pulsar  a 


- variable,  depends  on  the  day 


The  peak  flux  density  (which  is  also  the  minimum 
detectable  flux  density)  for  the  different  days  is  shown 
in  Table  3-5;  this  table  also  contains  the  values  of  the 
mean  pulse  energy  and  mean  flux  density  computed  using 
equations  (3-4),  (3-5)  and  the  values  for  the  pulse  width 
given  in  Table  3-6. 

No  known  measurements  of  flux  density  and  energy  at 
low  radiofrequencies  exist  for  pulsar  PSR204S-16.  The 
lowest  frequency  data  on  this  pulsar  are  the  80  MHz  obser- 
vations made  by  Slee  and  Hill  (9) . They  attempted  to 
observe  the  pulsar  on  3 nights,  failing  to  detect  it  and 
giving  an  upper  limit  for  the  pulse  energy  of 
<0.17*10  26  Jm  2Hz  1.  The  failure  to  detect  this  pulsar 
is  not  surprising  since  it  has  been  reported  to  have  large 


Is 
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intensity 


i,  being  often  undetectable  (28, 


Our  experiences  at  45  MHz  confirm  this  behavior  also  since 
it  was  possible  to  detect  this  pulsar  only  in  the  first 
nights  of  the  observing  period;  it  faded  out  for  the  rest 
of  the  nights. 

An  unconfirmed  pulse  energy  of  tl.9  Jm-^Hz~*  has  been 
given  by  Backer  at  34  MHz  (30).  The  lowest  frequency  at 
which  published  data  on  pulse  energy  is  available  is  410 
MHz;  the  energy  is  0.120*10”^**  Jm~^Hz~*  (29,  31).  Our 
measurement  at  45  MHz  seems  rather  high  (perhaps  suggesting 
an  undetected  systematic  error) , although  no  direct 
comparison  is  possible  due  to  the  lack  of  published  informa- 
tion at  low  frequencies  on  this  pulsar. 

Pulse  energy  as  a function  of  the  frequency  has  been 
plotted  in  Figure  3-32,  including  the  data  mentioned  in 
the  last  paragraphs. 

3.2.b  Pulse  Width  for  Pulsar  PSR2045-16 

In  order  to  obtain  the  true  half  intensity  pulse  width 
of  PSR2045-16  at  45  MHz,  the  measured  half  intensity  pulse 
width  has  been  corrected  for  pulse  dispersion  and  post 
detection  time  constant  using  equations  (3-8)  and  (3-9) 
respectively. 

In  order  to  correct  for  pulse  dispersion,  the  follow- 
ing values  for  the  parameters  were  used: 


= 11.5  pc 


The  correction  for  pulse  dispersion 


the  above  values,  is: 


The  correction  due  to  the  post  detection  tine  constant 
of  ic  = 11.25  ras  is: 


The  true  (or  corrected)  half  intensity  pulse  width 
Wq  ^ is  as  before  (equation  3-10) : 


The  result  for  the  4 days  is  shown  in  Table  3-6. 


Observed  and  corrected  half  intensity 
pulse  width  and  pulse  width  at  base 
(assuming  triangular  shape). 


May  15,  1979 
May  16,  1979 
May  18,  1979 
May  19,  1979 


Width  at  base 
(ms) 
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The  lowest  frequencies  at  which  pulse  shape  and  width 
measurements  have  been  given  for  this  pulsar  are  150  and 
151  MHz  (10,  22).  Lyne  et  al.  (22)  have  given  the  separa- 
tion between  subpulses  at  a number  of  frequencies;  the 
values  are  plotted  in  Figure  3-33.  The  separation  between 
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the  leading  and  central  subpulses  and  the  trailing  and 
central  subpulses  is  plotted  separately. 

As  in  a previously  mentioned  case,  the  trailing  pulse 
seems  to  decrease  very  rapidly  in  intensity  toward  lower 
frequencies  (22),  and  it  is  likely  that  at  45  MHz  it  no 
longer  exists  or  at  least  it  can  be  very  weak.  The  mean 
value  at  45  MHz  for  the  pulse  width  at  its  base  (obtained 
by  measuring  the  width  at  half  intensity  and  assuming  a 
triangular  shape)  of  61  ms  has  also  been  plotted  in  Figure 
3-33. 

The  plotted  value  for  the  pulse  width  at  its  base,  lies 
close  to  the  extrapolation  of  a straight  line  passing 
through  the  points  representing  the  separation  between  the 
leading  and  central  subpulses.  A possible  conclusion  of 
this  is  that  the  pulse  observed  at  45  MHz  is  composed  of 
two  unresolved  pulses  which  might  correspond  to  the  leading 
and  central  subpulses  at  higher  frequencies;  the  pulse 
width  at  its  base  would  correspond  to  the  separation  between 
the  two  subpulses. 

No  polarization  measurements  have  been  made  for  this 
pulsar  at  low  frequency.  The  lowest  frequency  at  which 
polarization  measurements  are  available  is  409  MHz  (32,  33) . 
The  mean  percentage  of  linear  polarization  at  409  MHz  is 
40%,  and  it  changes  during  the  course  of  the  pulse,  reaching 
about  70%  for  the  first  subpulse.  The  stability  of  the 
linear  polarization  position  angle  is  low  for  the  wings  of 
the  integrated  pulse  profile,  due  to  the  occasional 
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occurrence  of  pulses  with  orthogonal  polarization  with 
respect  to  the  mean,  but  is  high  at  the  center  of  the 
pulse  (32).  If  the  degree  of  linear  polarization  of  the 
pulse  is  higher  at  lower  frequencies,  there  will  probably 
be  a distortion  in  the  pulse  shape  when  it  is  received 
with  a linearly  polarized  antenna.  This  distortion  will 
depend  on  the  variation  of  the  polarization  angle  within 
the  pulse.  Assuming  that  the  linearly  polarized  component 
in  the  central  part  of  the  pulse  is  parallel  (or  nearly 
parallel)  to  the  antenna,  then  the  integrated  pulse  profile 
will  appear  strong  in  the  center;  however  the  slopes  of  the 
leading  and  trailing  edges  will  be  weaker  and  the  averaged 
pulse  profile  will  be  narrow.  More  details  about  the 
effect  of  the  degree  of  linearly  polarized  component  on  the 
pulse  profile  and  apparent  intensity  variations  are  given 
in  the  conclusions. 

3.2.C  Period  of  Pulsar  PSR2045-16 

from  Observations  on  Two  Consecutive  Days 

The  apparent  mean  pulse  period  of  a pulsar  can  be 
calculated  from  the  elapsed  time  between  two  pulses;  in 
order  to  increase  the  accuracy  it  is  desirable  that  the 
elapsed  time  between  the  two  pulses  be  much  longer  than 
the  pulsar  period.  An  attempt  was  made  to  obtain  the 
apparent  mean  pulse  period  from  observations  made  24*' 
apart  (elapsed  time  t24h)  for  pulsar  PSR2045-16.  On  two 
occasions  the  pulsar  was  detected  on  two  consecutive  days; 
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the  first  occasion  was  May  15  and  16,  1979,  and  the  second 
May  18  and  19,  1979. 

The  pulse  period  can  be  obtained  by  subtracting  the 
epoch  of  apparition  of  the  pulse  E on  two  consecutive  days 
and  dividing  by  an  integer  number  N. 

The  epoch  of  the  pulse  can  be  calculated  from  the 
relationship: 


Tstart  = t^me  start  digitizing 
n = bin  number  where  the  pulse  occurs 
R - plotting  time  resolution 

The  information  related  to  the  two  pairs  of  days  has 
been  summarized  in  Table  3-7. 

The  mean  apparent  period  can  be  calculated  from  the 
relationship: 


where  AE:  difference  in  epoch  of  pulse  apparition  on  two 

consecutive  days  (elapsed  time  between  pulse 
apparitions) 

The  mean  apparent  periods  obtained  with  this  method 


periods  in  Table 


May  15  and 


(3-17) 


Assuming  an  elapsed  time  of  24n  (86400  s) , a resolu- 
tion of  0.0315  s and  a period  of  1.961386  s,  the  following 
accuracy  for  the  period  is  obtained: 

Ac  = 0.36  us 

This  accuracy  should  bring  the  two  pulses  in  phase 
within  ±0.5  bin. 

It  is  not  clear  what  caused  the  discrepancies  between 
the  DOPVEL  period  and  the  period  deduced  from  the  observa- 
tions; it  could  be  an  instrumental  problem  or  a precision 
problem  in  the  calculations  in  the  DOPVEL  program.  Both 
possibilities  were  checked  but  the  result  was  inconclusive; 
further  investigations  need  to  be  made. 

3,2,d  Pulsar  PSR2Q45-16  Pulse  Intensity 

Figure  3-34  is  a plot  of  the  relative  pulse  intensity 
as  a function  of  time  for  the  pulse  obtained  on  May  16, 

1979;  this  pulse  is  the  strongest  and  best  defined.  Each 
point  plotted  represents  the  relative  intensity  with  respect 
to  the  noise  average,  averaged  over  10  periods;  the  time 
axis  scale  is  in  pulsar  periods.  It  can  be  seen  that  the 
pulse  has  large  intensity  variations.  Most  of  the  time  it 
has  a positive  value  with  respect  to  the  noise  average,  and 
it  seems  to  show  an  almost  periodic  intensity  modulation  of 
■vl40  periods  (peaks  at  periods  20,  160  and  320). 
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CHAPTER  4 
CONCLUSIONS 

Two  pulsars  have  been  successfully  detected  using 
narrow  bandwidth  and  long  averaging  times.  Pulsar 
PSR1133+16  was  detected  with  the  University  of  Florida  large 
array  at  26.3  MHz  and  pulsar  PSR2045-16  was  detected  using 
the  Maipu  Radio  Observatory  large  array  at  45  MHz.  Pulsar 
PSR1133+16  was  strong  enough  to  be  considered  statistically 
significant  (pulse  deflection  larger  than  3 standard  devia- 
tions) in  the  data  of  2 nights  (out  of  10  nights  reduced) . 
Pulses  of  PSR2045-16  were  found  on  4 nights  (out  of  10 
reduced)  with  enough  intensity  to  reach  the  3 standard  devia- 
tion level.  Typical  averaging  times  for  both  pulsars  were  in 
the  range  400  to  600  periods. 

Both  pulsars  show  large  intensity  variations  on  a time 
scale  of  days.  Typically,  one  would  appear  after  an  aver- 
aging time  as  short  as  •'-200  periods,  and  then  would  fade, 
being  undetectable  for  averaging  times  as  long  as  •'-1500 
periods  over  succeeding  days. 

The  peak  pulse  flux  density  obtained  for  PSR1133+16  is 
about  150  Jy;  for  PSR2045-16  the  peak  flux  density  is  in  the 
range  of  120  to  140  Jy.  These  values  (which  have  relatively 
large  probable  error,  in  the  order  of  50%),  are  higher  than 
the  values  of  the  peak  flux  density  measured  by  others  at 


nearby  frequencies.  They  are  also  larger  than  the 
extrapolated  flux  densities  from  published  data  at  higher 
frequencies. 

Also,  the  pulse  shape  and  half  intensity  width  seem 
to  differ  from  extrapolated  values  based  on  data  at  higher 
frequencies. 

In  order  to  find  a satisfactory  explanation  of  the 
intensity  variations  of  the  two  pulsars,  it  is  necessary 
to  consider  both  previously  reported  intensity  variations 
and  the  data  available  on  polarization  at  higher  frequencies. 

Pulsar  PSR1133+16  and  several  other  pulsars  have  been 
reported  to  show  extremely  variable  flux  densities  both  from 
pulse  to  pulse  and  on  a longer  time  scale,  when  observed  at 
81.5  MHz  (34).  These  intensity  variations  are  not  believed 
to  be  caused  by  interplanetary  scintillation,  and  therefore 
must  be  attributed  to  the  source. 

Pulsar  PSR2045-16  has  been  reported  as  often  being 
undetectable  for  several  days  at  408  MHz,  but  when  it  is 
active  it  can  produce  pulses  stronger  than  any  other  of  the 
first  pulsars  discovered  (28) . 

Besides  these  intensity  variations,  the  intensity  of 
the  received  signal  will  show  larger  and  more  frequent 
variations  when  received  by  a linearly  polarized  antenna, 
if  the  pulses  are  strongly  linearly  polarized.  Manchester 
et  al.  have  found  that  most  pulsars  present  a relatively 
high  degree  of  linear  polarization  (32) . PSR1133+16  was 

found  to  have  a degree  of  linear  polarization  of  33«  at 
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403  MHz,  50%  at  86  MHz  and  60%  at  63  MHz,  which  is  the 
lowest  frequency  at  which  polarization  has  been  measured 
for  this  pulsar  (35,  36).  There  is  a clear  dependence  of 
percentage  of  linear  polarization  upon  wavelength,  and  the 
tendency  is  for  the  percentage  of  linear  polarization  to 
increase  with  increasing  wavelength.  Pulsar  PSR1133+16 
has  a complex  pulse,  composed  of  two  subpulses;  the  per- 
centage of  linear  polarization  is  variable  throughout  the 
pulse  and  reaches  about  70%  in  the  middle  of  it  at  403  MHz. 
For  individual  pulses,  the  position  angle  of  the  plane  of 
linear  polarization  is  highly  stable  for  all  but  the  first 
component,  which  usually  shows  two  preferred  positions 
perpendicular  to  each  other.  The  result  of  this  is  that 
the  degree  of  linear  polarization  of  the  integrated  pulse  is 
low  for  the  first  subpulse  and  high  for  the  rest  (32) . 

Pulsar  PSR2045-16  has  a complex  pulse  composed  of  3 
subpulses.  Manchester  et  al.  (32)  measured  the  degree  of 
linear  polarization  to  be  40%  at  403  MHz,  with  the  polariza- 
tion being  even  larger  for  the  first  subpulse  (~70%) . The 
position  angle  of  the  plane  of  linear  polarization  varies 
smoothly  throughout  the  pulse;  a change  of  about  30°  occurs 
between  the  first  and  third  subpulses.  This  pattern  of 
position  angle  variation  within  a pulse  is  stable  for  many 
consecutive  pulses.  No  known  polarization  measurements  are 
available  at  lower  frequencies  for  this  pulsar. 

Even  if  the  linear  polarization  position  angle  were 
stable  when  measured  at  a point  above  the  terrestrial 
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ionosphere,  the  ionosphere  can  introduce  relatively  large 
variations  in  the  angle  as  measured  at  ground  level.  This 
is  due  to  variations  in  the  amount  of  Faraday  rotation  in 
the  ionosphere.  This  effect  is  particularly  strong  at 
longer  wavelengths  because  the  angle  of  rotation  0 is 
proportional  to  the  square  of  the  wavelength  X: 


where  RM  = rotation  measure  of  the  ionospheric  path. 

The  RM  depends  linearly  on  the  electron  density,  N, 
the  component  of  the  earth's  magnetic  field  intensity  which 
is  parallel  to  the  direction  of  propagation,  B„  , and  the 
thickness  of  the  ionosphere,  L: 


It  can  be  shown  that  the  percent  variation  of  0 is 
equal  to  the  percent  variation  of  the  electron  density  N, 
if  the  thickness  L and  the  magnetic  field  B,  of  the  ionosphere 
remain  constant;  that  is: 


It  can  also  be  shown  that  ^ equals  ^ when  N and  Bn 
are  constant,  and  equals  ^ when  N and  L are  constant. 

Manchester  has  obtained  values  of  ionospheric  RM 
ranging  from  0.1  to  6 rad  m depending  on  the  time  of  day 
and  the  declination  of  the  pulsars  (37) . Taking  these  two 
extreme  values,  it  is  possible  to  obtain  the  change  in 
rotation  angle,  A0,  for  let  us  say,  a 18  change  in  AN  at 


447° 


radiation  is  relatively  highly  linearly  polarized,  this 


should  produce  a significant  drop  in  the  intensity  of  the 
received  signal  when  using  a linearly  polarized  antenna. 
This  last  effect  can  provide  a possible  explanation  of  the 
disappearance  of  the  pulsar  signal  when  switching  to  a 
different  beam  of  the  26.3  MHz  array. 

Regarding  the  effect  of  the  polarization  throughout 
the  pulse,  the  pulse  shape  can  be  modified  when  using  a 
linearly  polarized  antenna;  the  signal  will  appear  more 
intense  for  the  part  of  the  pulse  with  polarization  that 
matches  the  polarization  of  the  antenna,  and  abnormally 
weak  for  the  part  of  the  pulse  with  position  angle  unfavor- 
able for  the  antenna.  Also  if  a given  part  of  the  pulse 
presents  random  changes  in  the  position  angle  from  normal 
to  orthogonal  (which  led  to  a low  integrated  polarization 
for  that  part  of  the  pulse) , the  intensity  of  the  averaged 
pulse  will  be  abnormally  weak  in  that  part.  This  seems  to 
be  the  case  with  PSR1133+16,  in  which  the  observed  intensity 
at  26.3  MHz  of  the  second  subpulse  is  usually  stronger  than 
the  intensity  of  the  first  subpulse;  the  latter  is  often 
undetectable  or  appears  only  occasionally.  The  separation 
of  these  two  subpulses  at  26.3  MHz  is  much  larger  (M25  ms) 
than  the  predicted  separation  for  the  two  subpulses  (which 
compose  the  main  pulse)  obtained  from  data  at  higher  fre- 
quencies; however  the  separation  at  26.3  MHz  does  agree 


results  obtained 
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The  values  for  the  peak  flux  density,  the  mean  flux 
density  and  the  mean  energy  are  higher  than  the  values 
obtained  from  extrapolation  at  higher  frequencies,  with 
the  exception  of  a value  of  peak  flux  density  of  150  Jy 
obtained  at  the  relatively  high  frequency  of  111.5  MHz  for 
PSR1133+16.  Izvekova  et  al.  (7)  and  Kuz'min  etal.  (20) 
made  simultaneous  observations  of  this  pulsar  at  several 
frequencies  in  order  to  obtain  the  spectrum.  They  found 
that  the  average  spectrum  (average  over  several  nights) 
presents  a turn-over  at  about  60  MHz.  Nevertheless  they 
also  found  that  the  spectrum  shows  large  variations  from 
day  to  day  in  the  range  61  - 151  MHz,  often  with  a change  in 
the  sign  of  the  spectral  index.  It  is  interesting  to  note 
that  the  spectrum  of  some  individual  nights  does  not  show  a 
turn-over  but  rather  an  increasing  intensity  with  decreasing 
frequencies.  This  might  explain  (at  least  in  part)  the 
unusually  large  intensity  of  the  pulses  at  26.3  MHz,  when 
compared  with  average  values  obtained  by  other  authors  at 
higher  frequencies. 

Averaging  of  pulses  for  two  different  nights  was  not 
possible.  For  pulsar  PSR2045-16,  the  period  obtained  from 
the  DOPVEL  program  differs  from  the  period  deduced  from 
observations  on  two  consecutive  nights  by  amounts  between 
8 to  28  us;  this  discrepancy  is  large  compared  with  the 
required  precision  for  averaging  over  two  consecutive  days 
(i.e.,  t tenths  of  us).  No  attempt  was  made  for  PSR1133+16 
to  average  data  of  two  different  nights  since  the  only  two 


are  27  days 


apart.  It  is  not  clear  up  to 
the  present  whether  the  discrepancy  of  the  periods  is  due 
to  instrumental  problems  or  insufficient  precision  in  the 
calculations  of  the  predicted  period. 

From  the  results  obtained  in  this  thesis,  it  has  been 
shown  that  it  is  possible  to  detect  pulsars  at  both  26.3 
MHz  and  45  MHz,  with  certain  limitations.  In  the  light  of 
the  results  obtained,  some  improvements  and  modifications 
in  the  equipment  and  data  reduction  can  be  suggested  which 
will  make  it  possible  to  obtain  even  more  information  about 
the  received  emission. 


One  such  modification  is  already  being  put  into  effect. 

A greatly  improved  transistorized  26.3  MHz  receiver  which 
will  replace  the  old  receiver  is  under  construction  and  is 
nearly  completed. 

Recording  at  a much  wider  bandwidth  (perhaps  50-100  kHz) , 
and  utilizing  a de-dispersion  technique,  would  considerably 
increase  the  signal-to-noise  ratio.  The  wide  band  signal 
would  initially  be  recorded  using  an  instrumental  tape 
recorder  of  suitable  bandwidth.  It  could  then  be  played  back 
into  the  Saicor  Real  Time  Spectrum  Analyzer,  the  output  of 
which  could  be  digitized  and  further  analyzed  by  means  of  a 
microprocessor  or  a minicomputer  (like  the  PDP  1134). 

After  the  data  has  been  converted  to  a power  spectrum  appear- 
ing in  the  200  output  frequency  channels  of  the  Saicor 
spectrum  analyzer,  the  frequency  dispersion  can  then  be 
removed  by  means  of  '-1--  - 


the  computer. 


s-bandwidth 
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system  of  this  type,  the  signal-to-noise  ratio  would  be 
increased  by  a factor  of  , or  3.5,  over  what  it  was 

with  the  8 kllz-bandwidth  system. 

Another  advantage  of  a wide  bandwidth  system,  and  per- 
haps an  even  more  important  one,  is  its  ability  to  deter- 
mine the  degree  of  linear  polarization  of  the  pulsar  signals 
by  making  use  of  Faraday  rotation  in  the  terrestrial  iono- 
sphere. As  stated  previously,  the  rotation  of  the  plane  of 
linear  polarization  varies,  due  both  to  ionospheric  fluctua- 
tions and  to  a continual  change  in  the  ray  path  as  the  earth 
rotates,  so  that  the  polarization  plane  is  frequently 
parallel  and  frequently  perpendicular  to  the  antenna  dipoles. 
The  degree  of  linear  polarization  is  readily  calculable  from 
the  measured  intensities  at  times  of  parallel  and  perpendicu- 
lar orientations.  This  appears  to  be  an  excellent  method  for 
determining  linear  polarization,-  it  has  apparently  not  pre- 
viously been  used  for  pulsars. 


Also,  if  the  rotation  measure  (RM)  of  the  ionosphere 
above  the  observatory  can  be  determined  independently,  the 
position  angle  of  the  plane  of  linear  polarization  of  the 
pulsar  signal  at  the  top  of  the  ionosphere  can  be  calculated. 

It  may  also  be  possible  to  obtain  still  more  informa- 
tion about  polarization.  Since  the  dipoles  of  the  26.3  MHz 
array  in  the  Dixie  County  radio  observatory  have  the 
capability  to  be  rotated  to  any  position,  half  of  the  array 
can  be  rotated  in  90°  respect  to  the  remaining  half,  making 
possible  to  obtain  information  about  the  circularly 


polarized  component  of  the  emission.  If  for  some 
there  is  enough  signal-to-noise  ratio,  it  would  also  be 
possible  to  set  the  antenna  in  such  a way  that  one  third 
of  the  array  is  oriented  at  0°,  one  third  at  45°  and  the  last 
third  at  90°;  from  the  information  obtained  in  this  configu- 
ration it  is  possible  to  deduce  the  Stokes  parameters  of  the 
incident  wave.  Any  pulsar  polarization  measurement  made  at 
26.3  MHz  and  45  MHz  would  be  useful  since  none  have  been  made 
at  frequencies  lower  than  63  MHz  (36). 

Improvements  need  to  be  made  in  the  calibrations  in  order 
to  obtain  more  aocurate  values  for  the  peak  flux  density. 

This  can  be  accomplished  by  simulating  the  pulsar  pulse  using 
a precise  time  standard  and  a frequency  synthesizer,  in  order 
to  obtain  periods  within  some  p seconds  of  the  observed  period. 
It  should  be  possible  to  obtain  estimations  of  peak  flux 
densities  by  digitizing  and  processing  this  signal  in  a similar 
manner  to  the  actual  pulsar.  It  will  be  also  necessary  to 
find  out  the  reason  for  the  discrepancy  between  the  DOPVEL 
period  and  the  period  obtained  from  the  observations;  simula- 
tions with  an  accurate  period  should  help  in  determining  if 
the  discrepancy  is  due  to  instrumental  problems.  A more  care- 
ful analysis  of  the  accuracy  in  the  calculations  of  the 
DOPVEL  program  needs  to  be  made.  The  solution  of  this  problem 
will  make  it  possible  to  average  several  nights  and  obtain  the 
long-term  average  pulse  characteristics. 


APPENDIX  A 
LIST  OF  PULSARS 

Tables  A-l  and  A-2  contain  the  pulsar  candidates  to 
observed  using  the  University  of  Florida  26.3  MHz  array 
the  45  MHz  array  at  the  Maipu  Radio  Observatory. 

Table  A- 3 gives  the  references  from  which  the  data 
the  different  frequencies  were  obtained. 
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Table 
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408 


List  of  references  from  which  the  data  at 
the  different  frequencies  was  obtained. 
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APPENDIX  B 

SAMPLING  RATE,  TIMING  PULSE  FREQUENCY 
AND  DATA  POINTS  PER  SECOND 

For  the  configuration  of  the  data  reduction  system  used 
(Figure  2-2),  the  sampled  values  are  sent  to  the  Amdahl  470/V6 
computer  in  real  time.  In  that  condition  the  rate  of  trans- 
mission of  data  into  the  Amdahl  is  the  limiting  factor  for  the 
sampling  rate. 

The  absolute  maximum  transmission  rate  to  the  Amdahl 
470/V6  is  1200  Band  (120  characters/second) . The  maximum  use- 
ful rate  is  60  characters/second;  40  characters/second  is  a 
more  conservative  (but  reliable)  rate.  Since  the  A/D  con- 
verter has  12  bits,  which  means  3 characters  per  data  point, 
the  transmission  rate  will  be  40/3  = 13  data  points  per 
second . 

The  number  of  data  points  per  second  (DPS)  obtained  from 
the  KIM  I is  given  by  the  relationship  (39): 

DPS  = TPP»NAVExTSAMP  ,B_1> 

where  TPP  = KIM  I timing  pulse  period  (ms) 

NAVE  = Number  of  sampled  values  averaged  before 
sending  to  the  terminal 

TSAMP  = Parameter  for  the  KIM  I program 

The  TPP  for  the  KIM  I must  not  be  shorter  than  4 ms 
(250  Hz) ; the  value  used  in  the  digitizing  process  for  the 
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pulsar  was  8 ms  (125  Hz) . The  number  of  sampled  values 
averaged  (NAVE)  before  sending  through  the  terminal  to  the 
Amdahl  is  fixed  at  8. 

Table  B-l  gives  the  relationship  between  the  TPP,  TSAMP 


a function  of  T 


31.25  15.625  10.4167  7.8125 

15.075  7.8125  5.2083  3.906 


s Points  Per  Second  ( 


In  order  to  meet  the  transmission  rate  requirements  of 
13  data  points  per  second  using  8 ms  for  the  TPP,  a safe  value 
of  TSAMP  = 02  must  be  used;  that  gives  a value  of  7.8125  DPS. 

If  the  tapes  are  played  back  at  the  same  speed  that 
they  are  recorded,  then  the  resolution  would  be  also  7.815  DPS, 
which  is  a poor  resolution.  However,  if  the  tapes  are  played 
back  with  a slow-down  of  4:1  the  time  resolution  will  be  in- 
creased by  a factor  of  4. 

Table  B-2  gives  the  number  of  DPS  obtained  in  real  time 
or  original  pulsar  time  (39). 


Table  B-2.  Number  of  data  points  per  second  in  real 
and  original  pulsar  time. 


Real  time  15.675  10.4167  7.8125  6.250  5.2083 

Original 

time  62.7  41.667  31.25  25.00  20.8332 


Since  NAVE  is  fixed  at  8,  that  means  that  8 sampled  values 
are  averaged  together  before  producing  one  data  point;  in  that 
way.  the  tape  is  being  sampled  8 times  faster  than  the  data 
point  rate.  As  an  example,  if  a TPP  = 8 ms  and  a TSAM  = 02 
are  being  used,  7.8125  DPS  are  obtained,  but  the  tape  is  being 
sampled  at  the  rate  of  62.5  samples/ second  in  real  time,  which 
in  original  pulsar  time  is  250  samples/second. 


APPENDIX  C 
SCHEMATICS 
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APPENDIX  D 
COMPUTER  PROGRAMS 
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APPENDIX  E 

26.3  MHz  ARRAY  CALIBRATIONS 

In  order  to  obtain  the  minimum  detectable  flux  density 
and  then  estimate  the  peak  flux  density  for  pulsar  PSR1133+16, 
several  calibrations  were  performed  with  the  array. 

The  standard  noise  source  was  provided  by  2 Sylvania 
S722  noise  diodes;  as  a secondary  noise  source  a Hewlett 
Packard  amplifier  was  used.  Several  radio  sources  of  known 
flux  density  and  the  galactic  background  temperatures  were 
also  used  in  the  calibrations. 

Figure  E-l  is  a block  diagram  of  the  26.3  MHz  array 
signal  transmission  and  calibration  system.  Usually  the 
calibration  signals  are  introduced  at  point  © (Figure 
E-l).  In  order  to  improve  the  calibrations,  the  signal  was 
later  introduced  at  point  © (Figure  E-l);  by  doing  this, 
any  uncertainties  in  the  gain  (or  losses)  of  the  10  pre- 
amplifiers and  the  2 N-S  Butler  matrices  are  avoided. 

Two  1 input-to-6  outputs  50-ohm  power  splitters  were 
used  to  split  the  calibration  signal  into  the  10  preampli- 
fier inputs.  To  properly  match  all  the  outputs,  two  power 
splitters  output  were  loaded  with  50-ohm  resistors. 

As  a result  of  introducing  the  calibration  signal  at  © 
and  then  at  © , a net  gain  of  :6  dB  was  determined  for  the 
10  preamplifiers  plus  the  2 N-S  Butler  matrices;  several 
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other  gains  and  losses  were  also  determined  for  the 
preamplifiers,  250  KHz  filters  and  aluminum  coaxial  lines; 
all  the  losses  and  gains  determined  are  shown  in  Figure  E-l. 

E.l  Losses  Between  the  Dipoles  and  Point  (5) 

The  procedure  to  obtain  the  losses  between  the  dipoles 
and  point  (5) , which  includes  the  10  E-w  Butler  matrices, 
the  hybrid  rings,  the  matching  networks  and  the  parallel 
transmission  lines,  is  not  as  straightforward  as  it  is  to 
obtain  gain  and  losses  for  the  rest  of  the  system.  An 
indirect  method  was  therefore  used  in  this  thesis.  First, 
the  recorder  deflection  due  to  the  galactic  background  was 
measured.  Then  a signal  from  the  noise  generator  was  intro- 
duced at  point  © (using  the  indicated  switch) , and  was 
adjusted  to  give  the  same  recorder  deflection  as  that  from 
the  galactic  noise.  From  the  ratio  of  the  noise  temperature 
of  the  introduced  signal  to  that  of  the  part  of  the  sky  from 
which  the  galactic  signal  was  obtained,  which  is  available 
in  the  literature,  the  losses  can  be  obtained.  A summary  of 
the  results  is  shown  in  Table  E-l. 

The  noise  temperature  for  the  noise  diodes  has  been 
calculated  using  the  following  relationship  (17) : 


where  Tn  = ambient  temperature  (=  290°K) 
e = electron  charge  (=  1.6* 10” 19  C) 

X = diode’s  plate  direct  current 
R = plate  load  resistance 
k = Plank  constant  (=  1.38*10-23  J "K| 


126 


I.JB 


Substituting  in  equation  (E-l)  the  values  for  the 
constants,  using  the  value  R = 50  ohm  and  expressing  I in 
milliamperes,  the  following  relationship  can  be  deduced: 

T = 290  + 289.85  I (E-2) 

Equation  E-2  gives  the  temperature  in  °K  and  it  is  valid 
for  R = 50  ohm. 

Galactic  background  temperatures  for  the  first  two 
points  in  the  sky  have  been  obtained  from  Turtle  et  al. (40) 
at  26.3  MHz.  The  value  for  the  third  point  was  obtained  from 
measurements  made  by  Blythe  (26)  at  38  MHz  and  corrected  to 
26.3  MHz  using  the  galactic  background  spectral  index  given 
by  Findlay  (27)  of  n = 2.65,  which  is  valid  for  the  range 
20-400  MHz. 

Zt  is  interesting  to  compare  the  losses  obtained  in 
Table  E-l  with  the  losses  obtained  by  summing  the  estimated 
dB  attenuations  along  the  path  between  the  dipoles  and  point 
© . The  list  of  estimated  or  measured  losses  of  the  indi- 
vidual components  is  given  in  Table  E-2  (41). 

The  average  of  the  measured  losses  given  in  Table  E-l  is 
7.2  dB  which  agrees  well  with  the  losses  computed  from 


"2.K" 1 
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AT  for  a radio  source  of  known  flux  density  S.  The 
relationship  that  gives  the  effective  area  is: 


AT  is  the  change  in  antenna  temperature  referred  to 
the  dipoles,  corrected  for  losses  e between  point  © and 
the  dipoles.  If  AT^  is  the  change  in  temperature  measured 
at  point  © , then  the  actual  change  in  antenna  tempera- 


ATa=  — — (E-4) 

where  e = transmission  line  efficiency,  which  in  this 
case  considers  the  losses  between  point  © 
and  the  dipoles. 

From  the  calibrations,  the  following  values  were 
obtained  for  3C433: 

S = 267  jy 
AT^  = 295-K 
Ta  = 25,000°I< 
e = 0.19 

Using  equations  E-3  and  E-4,  the  following  value  is 
obtained: 


This  value  needs  t 


e corrected  by  a factor  p.  = 0 


2 Array  program  (11) , which  considers  t 


The  corrected 
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The  value  of  the  effective  area  can  be  also  expressed 
in  terms  of  the  change  in  antenna  temperature  in  °K,  divided 
by  the  flux  density  in  Jy  (">K/Jy);  this  quantity  is  usually 


Substituting  in  E-5  the  values 
ATa  = AT^/e  = 1553«K 
S = 267  Jy 
Px  = 0.69 

the  following  value  for  A*  is  obtained 
A|  = 8.3  °K/Jy 
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